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ABSTRACT 

Kinetic energy stored in ISM bulk/ turbulent motions is a crucial ingredient to properly 
describe most properties of observed galaxies. By using Monte Carlo simulations, we 
investigate how this energy is injected by supernovae and dissipated via cloud collisions 
and derive the corresponding ISM velocity probability distribution function (PDF). 
The functional form of the PDF for the modulus of the velocity dispersion is 

p(v) oc v 2 exp[— {v/ay]. 

The power-law index of the PDF depends only on the value of the average cloud 
collision elasticity (e) as (3 = 2exp((e) — 1). If (3 and the gas velocity dispersion 
a are known, the specific kinetic energy dissipated by collisions is found to be 
oc a 2 ln(2//3)/(/3 — 0.947); in steady state, this is equal to the energy input from SNe. 
We predict that in a multiphase, low metallicity [Z w 5 x 1O~ 3 Z ) ISM the PDF 
should be close to a Maxwellian (J3 — 2) with velocity dispersion a ^ 11 km s _1 ; 
in more metal rich systems (Z £ 5 x 1O~ 2 Z ), instead, we expect to observe almost 
exponential PDFs. This is in good agreement with a number of observations that we 
review and might explain the different star formation modes seen in dwarfs and spiral 
galaxies. 
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1 INTRODUCTION 

An increasing number of recent studies have come to the 
conclusion that a large, if not dominant, fraction of the en- 
ergy injected into the interstellar medium (ISM) in various 
forms by stellar activity, resides in bulk/turbulent motions 
of the gas (e.g. McKee 1990, Lockman & Gehman 1991, Fer- 
rara 1993, Vazquez- Semadeni et al. 1995, Norman & Ferrara 
1996 [NF]). Energy in this form is required, for example, to 
reproduce the observed vertical distribution of the HI in the 
Galactic disk. The corresponding turbulent pressure is esti- 
mated to be from ^ 8 (McKee 1990) to ~ 30 (NF) times 
larger than the thermal one. Note that such large ratios 
imply a high porosity factor, Q, of the hot gas produced 
by supernova (SN) explosions: McKee (1990) estimates that 
Q ~ l.l(P turb /P tot ) i/3 . This dynamical property of the ISM 
has a dramatic importance in shaping the resulting struc- 
ture of the parent galaxy. A number of numerical models 
(Vazquez-Semadeni et al. 1995, 1998; Ballesteros-Paredes 
et al. 1999; Korpi et al. 1999) have been developed in or- 
der to describe the global properties of a turbulent ISM in 



the Galaxy improving on the pioneering works by Rosen 
et al. (1993, 1996) and Rosen & Bregman (1995). Although 
these studies provide a first approach to the study of the ISM 
properties on galactic scales, finite resolution limits their 
ability to properly describe shocks arising both as a conse- 
quence of SN explosions and cloud collisions, thus making 
their estimates of the kinetic energy dissipation rather un- 
certain. The same caveat should be made for the treatment 
of thermal instabilities (and hence cloud formation) which 
requires to spatially resolve the "Field length" (Hennebelle 
& Perault 1999). Galaxy formation studies in a cosmologi- 
cal context (Katz 1992, Mihos & Hernquist 1994, Navarro 
& Steinmetz 2000, Springel 2000), have by now firmly es- 
tablished that star formation cannot be properly regulated 
by forcing a phase transition of the gas into the hot phase 
by thermal energy input from SN explosions alone. Current 
computational power limits the spatial and the time scale 
resolution of the simulations, implying a too highly dissipa- 
tive treatment. Moreover small-scale physical processes usu- 
ally cannot be modeled with sufficient detail. A surprisingly 
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better description of galactic disk properties is obtained once 
an heuristic equation of state for the ISM is assumed with 
parameters adjusted to simulate the effects of a turbulent 
pressure. In these models a fraction of the available energy is 
invested in modifying the kinetic energy of the gas (Navarro 
& Steinmetz 2000). The recipe for this approach is relatively 
simple and it constitutes a promising path to understand 
galaxy formation, but it needs to be calibrated against well 
known situations as for example the Galactic ISM, where 
a large number of high quality observations on the HI dy- 
namics are available. Finally, Ferrara & Tolstoy (2000) point 
out that if energy injection is mainly regulated by massive 
stars, then the observed metallicity range of galaxies must 
be consistent with what inferred from the ISM kinetic energy 
budget. Along these lines they find that a tight relation be- 
tween gas velocity dispersion and metallicity in dwarfs must 
exist . 

In view of the above situation, we have decided to ex- 
plore how kinetic energy is injected by SN explosions and 
radiated away as converging gas flows collide at supersonic 
velocities. The present work is mainly building on previous 
investigations by our group which have concentrated on the 
analysis of high-res hydro- and MHD simulations of cloud 
collisions (Ricotti, Ferrara & Miniati 1997 [RFM], Miniati et 
al. 1997, 1999). The main effort there was devoted to quan- 
tify the amount of kinetic energy dissipation as a function of 
the collision parameters and cloud properties. Here we take 
a more global, although simplified, approach to understand 
the interplay between the two above mentioned processes 
(SN explosion and cloud-cloud collisions) regulating the dy- 
namics/energetics of the turbulent ISM. To this aim, we will 
adopt a statistical description based on a Monte Carlo ap- 
proach that allows us to derive a probability distribution 
function (PDF) for the ISM velocity dispersions. The PDF 
is then shown to hold directly measurable and easily inter- 
preted information about the kinetic energy regulation in 
the ISM. 

The paper is organized in the following manner. In § ^ 
we discuss our statistical approach to the modeling of a tur- 
bulent ISM and the model assumptions. In § ^ we present 
the results. In § ^| we show that our results are in agreement 
and can explain a number of observational results, that we 
review, from published literature. In § Owe justify the sim- 
ple assumptions of our model and in § g we summarize the 
results. 



2 MODEL ASSUMPTIONS 

We consider a galaxy of total HI mass Mc jv u distributed in 
clouds of mass M c and cloud filling factor f c . The clouds con- 
stituting the Galactic Cold Neutral Medium (CNM) phase 
have typically a gas number density n c — p c / p,m p ~ 50 
cm -3 , where p is the mean molecular weight and m p the 
proton mass. If the clouds are further idealized as homoge- 
neous spheres of radius R c , the total number of clouds is 
Ntot = Mcnm /M c . The cross section for cloud collisions is 
A — BitRc, where the unknown parameter B ~ 1 takes into 
account the effect of magnetic field, tidal forces and gravita- 
tional focusing, effective at low (5-6 km s _1 ) cloud relative 
velocities (Lockman & Gehman, 1991). It follows that the 
number density of clouds is Af c = 3A r tot / c /47ri?i?, and their 



collisional frequency v = J\f c Av r , where v r is the mean rel- 
ative velocity of colliding clouds. 

We use a set of Monte Carlo simulations to model the 
dynamics/dissipation of cloud collisions, and the energy in- 
put from SN explosions and stellar winds. We assume a ho- 
mogeneous and isotropic spatial distribution of the clouds 
and energy sources. In this case, the system is completely 
characterized by the evolution of the velocity probability 
distribution function (PDF) for the modulus of the velocity 
of the clouds, which in case of a Maxwellian distribution can 
be expressed as, 

p(v) oc v 2 exp (^-^j ' W 

where the 3D velocity dispersion, a, is related to the ob- 
served ID velocity dispersion by the relationship, a = 

Starting from an initial PDF (typically a Maxwellian), 
we evolve the distribution at subsequent time steps At — 
(vNtot)' 1 as follows. At each time step we select two clouds 
(we do not consider three body encounters) from the velocity 
distribution, with initial velocities V\ t i and V2,i, by randomly 
sampling the PDF and calculate the final velocities Vij and 
V2j after the collision. The expression for the new velocities 
are 

v l,f = v r + v cm — 2v r V cm COS If | (2) 
v 2,f = v r + v cm + 2V r V cm COS If , (3) 

where the relative velocity after the collision and the center 
of mass velocities are, respectively: 

2 £ r 2 2 
V r = + v 2.i - 

2vi,i«2,i(sin cos <p + cos 6 sin </>cos ip)] ; (4) 

v lm = T^M + v 2,i + 

2^1,^2,8 (sin. 6 cos (p -f- cos 9 sin (f> cos tp)\ ; (5) 

e(v r ,i, Rc, Z) is the elasticity of the collision, with v r ,i being 
t he rel ative velocity before the collision (see discussion in 
§ 2.1.2) and Z the gas metallicity. For the definition of the 
angles 8, tp, cp see the sketch in Fig. |l|. 

Every 7 = (i/gjvAt)" 1 collisions, where vsn is the SN 
rate in the galaxy, an energy E — kEs m is added to the sys- 
tem; Es n = 10 erg is the total energy of the SN explosion 
and k = 0.02 is the fraction^] converted in kinetic energy of 
the engulfed clouds (Thornton et al. 1998). As long as cloud 
coagulation, cloud destructions or formations are neglected 
in the model, the following simple argument suggests that 
the cloud collision rate, v, is constant. The scale height of 
the gaseous disk in galaxies is H oc v r (McKee 1990, Fer- 
rara 1993). Since in our model we assume constant gas mass 
and cloud mass (i.e. , N to t = const., at least in a statistical 
sense) we have M c ocv^ .It follows from the definition of v, 
given at the beginning of this paragraph, that v is constant. 
If, in addition, we assume constant SN explosion rate we 

* The adopted value of the coefficient n is perhaps controversial 
because based on ID hydro-simulations. In a more realistic model 
dynamical instabilities and turbulence could play a crucial role. 
We do not worry about this uncertainty because the qualitative 
results of this paper are unaffected by a different choice of K. 
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Figure 1. Sketch of the collision geometry and angle definitions. 

have 7 = const, since At = (i>Ntot)~ and Ntot is constant. 
To summarize, the model is characterized by four constant 
parameters which must be appropriately chosen to describe 
the galaxy under study: M c , R c , Z and 7. Note that in a 
self-consistent global model of the ISM these four parame- 
ters are not independent but correlated. Here we are not in- 
terested in a global modeling of the ISM that would require 
a different approach to the problem using high resolution 
hydro- simulations. Therefore we rely on observational and 
theoretical results found in literature to determine a range 
of realistic values for those parameters. 

2.1 Energy sources and sinks 

2.1.1 Turbulent pumping by supernova explosions 

The main source of the kinetic energy required to support 
the observed HI bulk motions in the Galaxy is provided by 
SN explosions and stellar winds (McKee 1990, NF). The 
latter authors calculated the source function for interstellar 
turbulence by using a detailed formalism to account for the 
interaction between SN shock waves and the clouds, which 
are thereby accelerated and put into motion. For brevity, 
we only repeat here the main results relevant to this study, 
deferring the interested reader to NF for further details. 

To simplify our treatment we only consider primary 
shocks, i.e. shocks associated with the SN explosion itself, 
as a source of kinetic energy for the clouds. Therefore we 
neglect secondary shocks, i.e. primary shocks reflected off 
the clouds; as seen from Fig. 1 of NF , their contribution is 
less than 1% both for SNe and winds. 

The distribution function Pi^'Qj,) for primary shocks 
due to the i-th source, where i = 1 (i = 2) refers to SNe 
(winds), can be written as 

P t (M) = to^&M-™ = (6) 
where M is the Mach number of the shock with respect 



to the intercloud gas supposed here to be at T — 8000 K 
and V is the volume affected by the energy input. The pre- 
vious expression contains two parameters (Roi and at) which 
completely define the self-similar stage of the shock radius 
evolution for the given source. Adopting the canonical val- 
ues of NF, we take i? ;/pc = (70,70) and a t = (9/7,9/2). 
The values for the radii of different turbulence sources (i.e. , 
isolated and clustered SNe, Winds, HII regions) are given in 
Table 1 of NF. In turn, the specific value of 70 pc is taken 
from Cioffi, McKee & Bertschinger (1988). The resultant 
primary shock distribution in the ISM is then 

i 

From this distribution we randomly extract a value of 
the shock Mach number (or velocity, v s ); the velocity of 
the engulfed cloud before interaction with the shock, v c o, is 
extracted from the actual PDF. The velocity acquired by the 
cloud in the direction of the shock propagation is (McKee 
& Cowie 1975, Miesch & Zweibel 1994) 

Ves = ( ATl) > < 8 > 

where A = 4 for an adiabatic shock and \ = 100 is the 
canonical value for the density contrast between cloud and 
intercloud medium^} Thus the final velocity of the cloud 
after the interaction with the shock is 

= («c0 + v ls + ZVcsVcO cos , (9) 

where £ is the angle between the directions of the shock 
propagation and the cloud motion before the interaction. 

2.1.2 Kinetic energy dissipation in the collisions 

Since the pioneering work of Spitzer, interstellar cloud colli- 
sions have always been considered as perfectly inelastic be- 
cause, for typical Galactic clouds, the cooling time of the 
shocked gas is shorter than the characteristic time scale of 
the collision. Nevertheless, this assumption might be rather 
crude for small clouds and/or primordial galaxies where the 
gas metallicity is low and therefore the cooling time is longer. 
In a previous paper (RFM) , we have investigated the depen- 
dence of the collision elasticity, e cc , defined as the ratio of 
the final to the initial kinetic energy of the clouds, on the ve- 
locity and mass ratio of the colliding clouds, metallicity and 
magnetic field strength. The problem has been studied via 
high-res numerical simulations. In the case of head-on colli- 
sions, RFM derived a handy analytical relationship that has 
been shown to correctly approximate the results of numeri- 
cal simulations: 

„ 1 ^cool V r /ir\\ 

Tfcc = 1 - e cc = — oc — — — ; (10) 

6n c R c n c R c A(v r ,Z) 

^cooi is the column density of the post-shock radiative re- 
gion, and A [erg cm 3 s _1 ] is the gas cooling function. For 

T X = i s an average value for the density contrast which 
is essentially derived by imposing pressure equilibrium between 
clouds and intercloud medium, as in the classical two-phase model 
scheme of the ISM. For small clouds this ratio can vary; however, 
the cloud velocity has only a mild, \~ 1//2 , dependence on its value. 
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clarity we summarize here the relevant results of RFM: (i) 
the kinetic energy dissipation in cloud collisions is minimum 
(i.e. the collision elasticity is maximum) for a cloud relative 
velocity v r ~ 30 km s _1 ; (ii) the above minimum value is 
proportional ZR\, where Z is the metallicity and R c is the 
cloud size: the larger is ZR 2 C the more dissipative (i.e. , in- 
elastic) the collision will be; (Hi) in general, we find that 
the energy dissipation decreases when the magnetic field 
strength and mass ratio of the clouds are increased and the 
metallicity is decreased, respectively. 

In principle, molecular hydrogen cooling could be im- 
portant, especially in a low metallicity gas. However, H2 in 
diffuse gas is very easily photodissociated by stellar UV flux. 
Therefore we neglect to a first approximation the effects of 
such molecule on the elasticity. 

In RFM we derived the elasticity for face-on collisions. 
We can estimate the mean elasticity of cloud collisions due 
to off-center encounters by assuming that the overlapping 
cloud area experiences a face-on collision whether in the rest 
of the cloud volume kinetic energy is conserved. With this 
simple hypothesis we obtain 



„ = l_ e= J£. [ 2Ra (1 _ JL) db= !hl 
1 2R C J Q V 2R J 2 



(11) 



where r/ cc is the energy dissipated for face-on collisions. With 
this definition, perfectly inelastic 3D collisions have a mean 
elasticity e ~ 0.5. In Fig. ^ we show e as a function of 
the cloud relative velocity, v r , for clouds with metallicities 
Z = 1,0.1 and 0.01 Zq. We have assumed cloud tempera- 
ture T — 48K, cloud radius R c = 0.5 pc and cloud density 
given by eq. |l^, derived form the assumption of pressure 
equilibrium of the multiphase ISM. The effect of magnetic 
field (neglected here) is both to increase the cross section of 
the collision and the elasticity of the face-on collision. 



2.2 Selection of model parameters 

As previously mentioned the parameters of our model can 
be fixed according to the derived observational constraints. 
For the Milky Way the values we adopt here are derived 
from data collected in Thronson & Shull (1990). The cloud 
gas number density is n c « 50 cm -3 and the cloud radius is 
R c = 0.1 — 6 pc; the cloud masses are therefore in the range 
M c ~ 0.0045 - 1000 M Q . A typical atomic cloud has R c = 
0.5 pc and mass M c ~ 0.5 M . Clouds with M c > 1000 M 
are self-gravitating and we assume that they cannot remain 
in the atomic phasej^]. The value of 7, expressing the number 
of cloud collisions between two subsequent SN explosions is 
in the range 1 — 1000. The elasticity of the collision depends 
on the gas metallicity Z, the cloud radius R c , and relative 
velocity v r , according to eq. Finally, we have adopted a 
SN rate of vsn = 0.04 yr _1 . Note that once a value for 7 is 
fixed, the SN rate determines the time step of the simulation 
but does not affect the properties of the PDF at equilibrium. 



t Even clouds with M c ^ 1000 Mq are fully molecular only in 
dense knots and should remain atomic in the outermost regions. 
But we neglect this complication. 
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Figure 2. Elasticity of cloud collisions, e, as a function of the 
cloud relative velocity, v r , for clouds with metallicities Z = 1, 0.1 
and 0.01 Zq. By definition perfectly elastic collisions do not ra- 
diate away their initial kinetic energy because the shock prop- 
agating inside the clouds remains always adiabatic. The peak 
of the cloud collision elasticity at v r ~ 30 km s — 1 corresponds 
to maximum column density of the post-shock radiative region, 
^cool v r /A. See text for more explanation. 



3 RESULTS 

3.1 A Test Case: Elastic Collisions 

We have first tested our code on the perfectly elastic cloud 
collisions case. We have chosen several different initial con- 
ditions for the PDF, and in all cases, after a number of colli- 
sions comparable to the total number of clouds in the box, a 
stationary Maxwellian PDF is approached, as expected. The 
number of clouds and the total energy are conserved (in this 
case we set the energy input to zero) within less than 0.1%. 



3.2 Constant Elasticity 

We now assume that the initial PDF is Maxwellian with 
velocity dispersion <taj = 10 km s" 1 for all cases discussed 
below. We run several simulations exploring the effects of 
different values of 7, R c , n c and elasticity. The elasticity 
range is ^ e ^ 1 and in each simulation is set to a constant 
value unrelated to the physical properties of the clouds. 

The typical number of collisions required to reach the 
stationary regime depends on the mean kinetic energy per 
cloud mass, (E), at the time when the stationary phase is 
reached: the lower is (E), the larger is the number of col- 
lisions required to reach the steady state. An example of a 
Monte Carlo simulation is shown Fig. []. 

The shape of the stationary PDF is well described by 
the parametric function 
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Figure 3. Evolution of the cloud velocity distribution, N(v) = 
Ntotp{v), from a Monte Carlo simulation with e = 0.7,7 = 
110, R c = 4.17 pc, M c = 40 M . The dashed line is the initial 
PDF, a Maxwcllian distribution with <xjf = 10 km s ; the solid 
lines show the PDF at t = 25 and 100 Myr. For this run, a sta- 
tionary PDF is reached after about 50 Myr. 



p(v) = 







r(3//9)<T 



exp 



(12) 



whose maximum occurs at v max — (2//3) 13 a. Here F(:r) is the 
gamma function. The distribution has total kinetic energy 
per unit mass 



(E) 



r(5//3) 2 



3 1.05a 2 



2/3- 0.947 



if 1< 8 < 2 



(13) 



For the Maxwellian (6 = 2), we have (E) = (3/2)<rf f and 
(E) = (3/40)cr 2 if /3 = 1. In the range 1 < 3 < 2 the function 
r(5//3)/r(3//3) « 1.58/(/3 - 0.947). 

We have performed several simulations to explore the 
effect of changes in the model parameters. If the mass of the 
clouds is not too large (<; 1000 Mg), the PDF is given by 
eq. ([[2]), with the index 3 which is a function of the mean 
elasticity of the collisions only and independent on the de- 
tails of the energy injection. However, in the case of massive 
clouds each collision dissipates an energy comparable to the 
energy injected by a typical SN explosion. Therefore, at the 
steady state, the number of cloud collisions per SN explosion 
is about one and the functional shape of the PDF depends on 
the SN energy input prescription. In Fig. ^ we show the tight 
correlation of 3 to the elasticity. The solid line, showing the 
best fit to the simulation points, has the simple expression 

8 = 2e (E_1) , < e < 1. (14) 

For the perfectly inelastic case, according to eq. ([n]) e = 0.5, 
implying 8 ~ 1.2. 

A second conclusion can be drawn from our results. The 
mean kinetic energy per unit mass, (E), of the equilibrium 
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Figure 4. Shape of the PDF as a function of the elasticity. The 
parameter /3 in eq. ( p^ ) is a function of the elasticity e only. The 
solid line shows the curve j3 = 2 exp (e — 1): for elastic collisions 
(e = 1) the PDF is Gaussian, for inelastic collisions (e = 0.5) the 
PDF is close to exponential (/3 = 1.2). 



distribution is simply related to the energy input per unit 
mass per cloud collision, Ei = kEsn /"/M c , through the dis- 
sipation parameter r\ = 1 — e by the following equation: 



E i = r ] {E) = (l-e)(E) 



r(3/B) a 



(15) 



where we have used the two above relations. Eqs. (p^)-([l5[) 
are useful to estimate ISM parameters, {E),e and Ei from 
the observed shape of the PDF, characterized by 8 and a. 
The typical value of M c is a function of metallicity and will 
be derived in the next Section. 



3.3 Physical Elasticity 

We now use the more realistic prescription, eq. (|ll]) for the 
elasticity of cloud-cloud collisions. The solid lines in Fig. ^ 
show the isocontours of constant 8 (left panel) and mean 
elasticity (e) (right panel) as a function of Z and R c from the 
Monte Carlo simulations. The cloud gas density is n c = 50 
cm -3 and the temperature T = 48 K, the typical values for 
a two-phase ISM with Z = Z Q (see RFM). 

The value of 3 is found to be related to elasticity sim- 
ilarly to eq. |l4|, where now (e) is the mean elasticity. The 
latter is defined by 



(e) = dv r p(ar,P, v r )e(v r ,Z, n c ,R c ) 
8 = 2e ((e> ~ 1) 0.5 < (e) < 1 



(16) 



where p(a,B,v r ) is the probability distribution function for 
the modulus of the relative velocity (rPDF). The values of 
3 and (e), derived solving the system of eqs. dli|), are shown 
in Fig. H (dashed lines). 
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Figure 5. Isocontours of j3 (left) and mean elasticity (e) (right) as a function of gas metallicity Z and cloud size R c . The solid lines 
refer to Monte Carlo simulation results; the dashed lines are the analytic solutions of eqs. |H| The cloud density is n c = 50 cm -3 , the 
temperature T = 48 K and a = 10 km s — 1 . 



The ID rPDF is given by 

PlD(v r ) = / d,V! / dv2PlD{vi)piD{V2)8[v r - (V! - V 2 )] (17) 



Here, 5 is the Dirac delta function and pm(v) oc 
exp[— (v/aio) 13 ] is the normalized ID PDF. It is easy to 
solve eq. ( |l7| ) for the particular cases (3 = 2 and ,3 = 1. For 
(3 — 2 the ID rPDF is a Gaussian with a r = \pla. In the 
(3 — 1 case, the rPDF is 



Pio(o"id,/3 = l,v r ) = 



crin 



- exp 



Vr 



(18) 



For arbitrary values of (3 eq. ( |17[ ) has to be solved nu- 
merically. Eq. (|l^) shows that the functional form of the 
ID rPDF and therefore the rPDF p(a,(3,v r ) changes de- 
pending on (3. Nevertheless, using Monte Carlo simula- 
tions we find that the parametric function p(o,(3,v r ) oc 



v 2 exp[ 



'(v r /a r ) l3r ] with 



(1 + /3 _3 )2/? and (3 r 



0(1 + 0.24/3" d ) is a good fit to rPDF. We have used this 



parameterization to solve the system of eqs. (|16 

The mean elasticity, and therefore f3, depend on a. In 
Fig. we show (e), calculated from eq. (M), as a function 
of a for Z = 5 x 10~ 2 Z e and Z = 5 x 10~ J Z©. The trian- 
gles show the analogous relation obtained from Monte Carlo 

simulations. 

As discussed in § 2.1.2, cloud-cloud collisions are more 
elastic when their relative velocity is about 30 Km s -1 . This 
can be understood if we recall that the elasticity of the colli- 
sion is proportional to the column density of the post-shock 
radiative region, T, coo i oc v s t coo i oc wf/A(T), where v a is the 
shock velocity, t coo i the cooling time of the shocked gas and 
A(T) is the gas cooling function. For small collision veloc- 
ities the post-shock gas is cold (T < 10 4 K) therefore the 
cooling time is relatively long. In this regime the elasticity 
increases with the cloud relative velocity about as vj? . When 
the cloud relative velocity is high enough that the post-shock 



gas is heated above T ~ 2 x 10 4 K, the cooling time becomes 
very short because of the Lya cooling. Therefore the elastic- 
ity decreases steeply. The maximum of the mean elasticity 
(e) at a w 11 Km sT x is the consequence of the dependence 
of the cloud-cloud collision elasticity on the velocity of the 
clouds. In other words E coo ; is maximum when a m 11 Km 
s" 1 . 

The dependence of (e) on a indicates that if a < a c ~ 11 
km s _1 the PDF cannot reach a stationary state, i.e. an 
instability occurs. Physically, this is due to the fact that in 
this low-cr regime an increase (decrease) of a produces an in- 
crease (decrease) of (e) . As a result less (more) energy is dis- 
sipated and a increases (decreases) further. The simulation 
results shown for this a range do not reflect equilibrium val- 
ues but are given at a time t ~ 100 Myr; eventually, depend- 
ing on the parameters of the simulation the PDF relaxes ei- 
ther to the stable high-a branch or to a ~ 0. The time scale 
for this process is essentially given by the growth rate of the 
instability, which is particularly fast for a low metallicity gas 
and vanishes as Z is increased above Z w O.IZq . The growth 
rate of the instability is r _1 ~ \Ei/(E) — (rj)\v oc (1— (e))v, 
where we remind that v is the cloud collision rate. The tri- 
angles in the unstable region are evolved from the same 
Maxwellian PDF with om — 10 km s _1 but with differ- 
ent values for Ei. The fact that the triangles in the unstable 
region lay very close to the solid line implies that eq. ( |l5| ) 
is approximatively valid also when the system is not in a 
steady state. If E t < Ef lt ~ 121[r(5//3*)/r(3//3*] m(2//3*) 
km 2 s" 2 , where 3* = f3(a = 11) is a function of Z,M C 
and n c (see Fig. Q), the PDF will be unstable. Therefore, 
from Fig. ^ we conclude that in a multiphase, low metallic- 
ity (Z w 5 x 1O -3 Z ) ISM the PDF should be close to a 
Maxwellian (f3 = 2) with velocity dispersion a £ 11 km s _1 ; 
in more metal rich systems (Z XL 5 x 1Q~ 2 Zq) instead we ex- 
pect to observe an almost exponential tail of the PDF with 
(3 w 1.2. If Z £ 0.1, the instability of the PDF is not fast 
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enough to put strong constraints on the value of a. We will 
discuss in § ^ the observational implications of these results. 

Finally, we estimate the typical cloud mass as a function 
of the metallicity. As discussed in RFM, a multiphase ISM in 
thermal equilibrium can exist in a narrow range of thermal 
pressure. The pressure range depends on the gas metallic- 
ity but the temperature of the clouds (CNM) is T ~ 48 K, 
independent of the metallicity. This implies that exist a re- 
lationship between the cloud density n c and Z. If we assume 
a cloud temperature T = 48 K and pressure equilibrium of 
the multiphase ISM we find 



,(Z) ~ n c {Z = 1) exp[0.0834(Z 



-1/2 



!)]• 



(19) 



Adopting this relation we explore two different models: 

Model A Since diffuse clouds are not self gravitating, 
we assume that their mass is a small (constant) fraction of 
the Jeans mass Mj = (HOM )T 3/2 n c (Z)" 1/2 . The tem- 
perature of the cloud is approximately constant (T — 48 
K), therefore the mass is M c oc n 3 (Z)~ 1 ^ 2 and the radius 
R c oc n c (Z)~ 1/2 . 

Model B Alternatively, we can derive the cloud size by 
imposing that the Kb abundance in the cloud is negligible. 
If the cloud is optically thin to the H2 photodissociating 
radiation in the Lyman-Werner bands (i.e. , has a column 
density E c = n{Z)R c (5x 10 14 cm -2 )), the H2 abundance 
will be xh 2 ~ 10~ 9 if the UV background in the galaxy has 
values comparable to the Milky- Way. We therefore assume 
that the column density of the clouds is a (constant) fraction 
of 5 x 10 14 cm -2 . In this case we find the scaling laws: R c oc 
n c {Zy x and M c oc n c {Z)~ 2 . 

In Fig. [?] we show the maximum (at a = 11 km s _1 ) 
mean elasticity (e) as a function of the gas metallicity, Z. 
The three solid lines show the mean elasticity assuming that 
cloud masses scale with metallicity as M C (Z) = M C (Z — 
l)[n c {Z)/n c {Z = 1)]" 1/2 for M C (Z = 1) = 0.1,1,10 M 
(i.e. , Model A). The dashed lines show the analogous re- 
lation assuming M C (Z) = M C (Z = l)[n c (Z)/n c (Z = l)] -2 
for M C (Z = 1) = 0.1, 1, 10 M (i.e. , Model B). The mean 
elasticity increases when the gas metallicity decreases as a 
consequence of the lower cooling rate of metal-deficient gas. 
When the mass of the clouds is small the mean elasticity is 
larger because of the smaller cloud column density, E c (we 
have seen in § 2.1.2 that the cloud-cloud collision elastic- 
ity is inversely proportional to E c ). The difference between 
models A and B is also caused by the dependence of the 
cloud column density on the gas metallicity. In Model A 
E C (Z) oc n c (Z) 1/2 and in Model B £ C (Z) is constant. 



4 OBSERVATIONAL LINKS 

The results presented above allow to make a number of pre- 
dictions concerning the ISM velocity PDF which can be 
tested against observations. The first attempts to determine 
the PDF of the ISM date back to the early 50s. These stud- 
ies used optical absorption lines toward OB stars or HI 21 
cm emission or absorption line to measure the velocity PDF 
of local gas at distances d ^ 1 kpc. More recently the ques- 
tion of the velocity PDF has received much less attention. 
In the following, we will briefly summarize the available ob- 
servational results; we particularly emphasize the point that 
the large majority of them indicate that a exponential PDF 



A 
V 




10 20 

a[km s _1 ] 



Figure 6. Mean elasticity (e), as a function of cloud velocity 
dispersion a. The solid lines show the analytical result obtained 
by solving eq. (|l|) for Z = 5 X 1CT 2 Z and Z = 5 X 1(T 3 Z 
(small integration errors are responsible for the wavy appearance 
of the lower line). The triangles show the analogous relation as 
obtained from Monte Carlo simulations. It can be shown that if 
a < CT C ~ 11 km s" 1 , the PDF is unstable. 



at large a is a much better fit than a Gaussian one to the 
velocity data. 

Adams (1949) compiled a catalog including 230 Ca K 
absorption line velocities produced by "clouds" along the 
line of sight to nearby OB stars. He found 150 components 
towards 120 stars with d < 500 pc and 80 components to- 
wards 43 stars with d > 500 pc. Using such catalog, Blaauw 
(1952) and Takakubo (1967), showed that a single exponen- 
tial could fit the velocity PDF of local interstellar clouds bet- 
ter than a single Gaussian. Analyzing the same data, Huang 
(1950) and Kaplan (1954) preferred a fit with a d -1 func- 
tion; Siluk & Silk (1974) fitted the high-velocity tail of the 
PDF with a v~' A function. Another set of optical interstellar 
lines, toward 132 stars was observed by Munch (1957). He 
used a different method, the "doublet ratio method" (Ca 
H and K or NaD2/Dl). His conclusions confirm Blaauw's 
results. 

Mast & Goldstein(1970) observed the radial velocities 
of a sample of 268 HI 21 cm emission clouds. In order to 
clearly resolve the clouds in velocity, they chose a sample 
of high latitude line of sights and estimate the peculiar ve- 
locity distribution after correction for solar motion and dif- 
ferential galactic rotation. It is evident from their Fig. 4 
that an exponential fits better the data than a Gaussian. 
Falgarone and Lequeux (1973), using data from Radhakr- 
ishnan et al. (1972 a,b,c,d), Goss & Radhakrishnan (1972) 
and Hughes, Thompson & Colvin (1971) analyzed the PDF 
of the ISM in the solar neighborhood. Their results are in 
agreement with those of Mast & Goldstein(1970) but their 
sample is too limited to discriminate the detailed shape of 
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z/z 

Figure 7. Mean elasticity (e), as a function of the metallicity Z 
(normalized to the solar metallicity), for a = 11 km s — 1 . We as- 
sume that the clouds are in pressure equilibrium in a multiphase 
medium. This assumption imply that the clouds have temperature 
T = 48 K and density n c (Z) ~ n c {Z = 1) exp[0.0834(Z~ 1 /2 _ i)] 
{i.e. , P(Z) = n c T ~ P(Z = 1) exp[0.0834(Z" 1 / 2 - 1)], where 
P(Z = 1) = 2400 cm -3 K). Here we assume n c (Z = 1) = 
50 cm -3 . The solid lines show the mean elasticity assuming 
that cloud masses scale with metallicity as M C (Z) = M C (Z = 
l)[n c {Z)/n c (Z = 1)]~ 1/2 for M C (Z = 1) = 0.1,1,10 M Q 
(Model A). The dashed lines show the analogous relation assum- 
ing M C (Z) = M C (Z = l)[n c (Z)/n c {Z = 1)]~ 2 for M C (Z = 1) = 
0.1,1,10 M Q (Model B). The rationale for models A and B is 
explained in the text. 



the distribution. Crovisier (1978) studied the kinematics of 
nearby HI clouds using the Nancay 21-cm absorption survey 
(Crovisier el al. 1978). The sample consisted of about 300 
cloud radial velocities. They also confirmed that the ID PDF 
is close to exponential, although their conclusion suffers of 
possible contamination from a population of intermediate- 
velocity clouds. Dickey, Terzian, & Salpeter (1978) obtained 
the radial velocity distribution of clouds from the Arecibo 
21-cm emission/absorption survey. Unfortunately, the statis- 
tics is too small to clearly discriminate the functional form 
of the PDF. Kim et al. (1998) measured the PDF for mo- 
tions out of the plane of the Large Magellanic Cloud. The 
peculiar velocity PDF, on scales larger than 15 pc, is well 
fitted by an exponential. 

In a slightly different context, Miesch, Scalo & Bally 
(1999) estimated the PDF of molecular line centroid velocity 
fluctuations for several nearby molecular clouds with active 
internal star formation. The data consist of over 75,000 13 CO 
line profiles divided among 12 spatially and/or kinematically 
distinct regions. These regions range in size from less than 
1 to more than 40 pc and show substantially supersonic 
motions. They find that 3 regions (all in Mon R2) exhibit 



nearly Gaussian centroid PDFs, but the other 9 show nearly 
exponential PDFs. 

In spite of this impressive amount of evidences, many 
authors have continued to make the assumption when mod- 
eling or interpreting ISM data that the velocity PDF for 
the diffuse and molecular clouds is Gaussian. This bias is 
probably due to the expectations from incompressible tur- 
bulence models of the ISM, which predict that the PDF 
is nearly Gaussian. For systems whose metallicity is larger 
than about 1/20 of solar, our model predicts instead that 
the velocity PDF should approximate an exponential distri- 
bution. The physical reason for this discrepancy resides in 
the dissipative nature of cloud collisions. From our results, 
an exponential PDF can only be obtained if the ensemble 
averaged elasticity is 0.5, indicating perfectly inelastic colli- 
sions. Hence, kinetic energy stored in turbulent motions is 
efficiently converted in thermal energy and radiated away in 
shocks arising in colliding flows. The (supersonic) turbulence 
decays very rapidly and a corresponding high rate of energy 
deposition is required to stir the gas. This is in agreement 
with the findings of up-to-date numerical simulations in the 
context of star formation (e.g. Smith, MacLow & Heitsch 
2000; Padoan & Nordlund 1999). 

Finally, we would like to comment on the dynamical 
instability discussed above (see Fig. ^) . If values of a corre- 
sponding to the unstable a < 11 km s _1 region are deduced 
from observations of a low metallicity system, because of 
the runaway character of the instability, a should rapidly 
decrease or increase assuming a constant SN energy injec- 
tion rate. If a decreases, leading to a loss of pressure support 
and hence to overall contraction, is likely to drive a similarly 
rapid increase of the star formation rate and energy injec- 
tion from SN explosions. This will cause a to start increas- 
ing. An increasing a will lead to overall expansion of the 
ISM and the star formation rate and the SN energy injec- 
tion will decrease back to the the starting value. Therefore 
a starts decreasing again unless the energy injection rate is 
sufficient to maintain a > 11 km s _1 . As a result a "bursty" 
star formation mode will occur, unless a is increased up to 
the stable region, thus stabilizing the star formation rate 
on a higher level. The period of the burst cycle should in- 
crease with the gas metallicity because the timescale of the 
instability does. A "bursty" mode of star formation has in- 
deed been observed in some low surface brightness (LSB) 
dwarf galaxies (Schombert, McGaugh & Eder 2001). Thus, 
the dynamical instability occurring together with the dissi- 
pation of kinetic energy investigated here might provide an 
interesting self-regulation of the star formation activity in 
low-metallicity and low surface brightness galaxies. The sug- 
gested self-regulation mechanism is rather speculative given 
the simplicity of our model and should be tested using more 
realistic numerical simulations. 



5 DISCUSSION 

Turbulent supersonic flows and gas self-gravitation create 
structures in the ISM that are very complex and hierar- 
chically structured. It is still under debate whether inter- 
stellar clouds are only transient phenomena produced by 
dynamical fluctuations or can be represented as long-living 
clouds. Numerical simulations (Wada, 2001) and observa- 
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tions (Schwarz & van Woerden, 1974) show that the mor- 
phology of the atomic component of the ISM is a complex 
network of filamentary structures. The filaments are proba- 
bly produced by oblique shock collisions and shear from dif- 
ferential galactic rotation. The typical size of the filaments 
in the simulations is about 1 pc, consistent with the ob- 
served typical size of atomic interstellar clouds. Contrary to 
self-gravitating molecular clouds, the atomic clouds do not 
show evident substructures. 

In our simple toy model clouds do not coalesce or are 
destroyed, and their mass, size and number density remain 
constant. Our assumptions are not physically justified but 
are the correct choice for a statistical description of the ISM. 

The outcome of slow inelastic collisions could produce 
cloud coalescence. If this process is important, eventually the 
cloud mass reaches the Jeans mass and the atomic cloud be- 
comes a self-gravitating molecular cloud. Perhaps this pro- 
cess, in conjunction with fragmentation of dense gas shells 
produces by SN explosions, is responsible for the regulation 
of the gas ratio between the atomic and molecular phases 
of the ISM. Simulations show that small clouds are prob- 
ably destroyed if they are passed by supersonic shocks. In 
this process mass is exchanged between atomic clouds and 
the intercloud diffuse gas. Indeed from numerical simula- 
tions it is evident that clouds are not long lived entities but 
instead they are continuously formed and destroyed. But 
despite the complex and time dependent morphology of the 
ISM, statistical analysis of the cloud properties in numeri- 
cal simulations shows that the typical cloud size, cloud mass, 
and cloud number density have constant equilibrium values 
(Wada & Norman, 2001). Therefore our model is not phys- 
ically realistic but our simple assumptions are justified in a 
statistical sense. For this reason we cannot study the ISM 
morphology or the physics that regulates mass exchange be- 
tween different phases of the ISM. But we believe that our 
model is reliable for a statistical study of kinetic energy dis- 
sipation in "cloud" collisions. At least on scales where the 
isotropy assumption in the ISM holds. 

Because of the aforementioned arguments, we suspect 
that adding complexity to the model (for instance assuming 
non homogeneous density clouds or a cloud mass spectrum) 
would not help substantially in making the results more 
realistic and reliable. Nevertheless, in a work currently in 
preparation, we are using 2D high-resolution hydrodynamic 
simulations to model the ISM in galaxies with primordial 
(metal-poor) composition. These simulations should be able 
to test the quality of the results presented in this paper. 

Finally we discuss to which extent the results of this 
work can be extended to different systems, namely sub- 
clumps in molecular clouds and when the main turbulent 
energy source is different from SN explosions. 

Observations of the velocity PDF in molecular cloud 
cores (Miesch, Scalo & Bally, 1999) have shown that the ID 
PDF is better represented by an exponential rather than 
Gaussian distribution. Based on the results of our model, 
that applies to atomic clouds, we have suggested that ki- 
netic energy dissipation in collision of sub-clumps is at the 
origin of the exponential velocity distribution. Even if we 
believe that an analogy between atomic clouds in the ISM 
and sub clumps in molecular clouds is possible, we have to 
keep in mind important differences. Molecular clouds, are 
self gravitating and harbor smaller clumps for which size- 



mass-velocity relations have been derived from observations 
(e.g., Larson, 1981). The inter-clump gas can be atomic es- 
pecially in the outermost regions. Our assumption of homo- 
geneity and isothermality for clouds of same sizes and the de- 
rived scalings can serve as a first-order approach to the real 
hierarchically structured molecular cloud. More important, 
in our model we do not include molecular chemistry and 
cooling. Note that, if molecular hydrogen is not dissociated, 
clump collisions are inelastic (because of radiative cooling 
from collisionally excited H2 rotational and vibrational lev- 
els) even if the gas has primordial composition (metal-free) . 

One of the results discussed in § H is that the shape of 
the PDF does not depend on the kinetic energy input pre- 
scription if the energy input event (SN explosion in our case) 
is much larger than the energy dissipated by a typical colli- 
sion event. This means that our results should remain valid 
even if the main kinetic energy source that fuels the ISM 
turbulence is not SN explosions. Stellar winds are are im- 
portant energy sources. Chimney and superbubbles created 
by multiple SN explosions in a stellar cluster can reduce the 
amount of energy input in the ISM and produce, instead, gas 
outflows from the ISM to the intergalactic medium. Finally 
Wada & Norman (1999) have shown that the rotational en- 
ergy of the disk is converted into turbulent energy of the 
ISM because of the combined effects of gas self-gravity and 
disk differential rotation. 



6 SUMMARY 

We have investigated the kinetic energy of the ISM in galax- 
ies as a function of their gas metallicity. We use a simple 
closed model where the energy injection from SN explosions 
balance the energy dissipated in cloud collision. Using Monte 
Carlo simulations, coupled with a simple (but statistically 
motivated) model for the clouds and previous numerical re- 
sults on the elasticity of collisions derived by RFM, we have 
been able to link the properties of the velocity PDF with the 
ISM kinetic energy dissipation. In particular, we find that: 

• The slope, f3, of the velocity PDF only depends on 
the value of the ensemble averaged elasticity (e) as f3 — 
2exp«e>-l). 

• The knowledge of f3 and of the velocity dispersion 
a of the gas allows to determine the specific kinetic en- 
ergy input from SNe, Ei through the simple relation Ei ta 
1.58cr 2 ln(2//3)/(/3 - 0.947); in steady state, this input is 
equal to the energy dissipated by collisions. 

• We predict that (see Fig. ^ in a multiphase, low 
metallicity (Z « 5 x W~ 3 Z Q ) ISM the PDF should be close 
to a Maxwellian (/3 = 2) with velocity dispersion a £ 11 
km s _1 ; in more metal rich systems (Z £ 5 x W~ 2 Zq), in- 
stead we expect to observe exponential PDFs with /3 « 1.2. 

• We have pointed out that, contrary to what is usually 
assumed, the available data firmly allow to conclude that the 
Galactic ISM velocity PDF is not Gaussian but has rather 
an exponential shape. This is in perfect agreement with our 
predictions. 

• Low metallicity systems with a < 11 km s" 1 are dy- 
namically unstable on time scales r oc [(1 — (e))v]~ l where v 
is the cloud collision rate. This systems are probably under- 
going a "bursty" mode of star formation. Indeed, some ob- 
servations suggest that LSB dwarf galaxies show signs of re- 
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peated and weak bursts of star formation which shape their 
irregular morphology appearance (Schombert, McGaugh & 
Eder 2001). On the contrary, a quiescent star formation 
mode as inferred for the Milky Way, should be associated 
with an ISM velocity dispersion a w 11 km s" 1 and a rela- 
tively high star formation rate. 
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